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Abstract 

b : 

. . . , We calculate both connected and disconnected contribution to the ir-N 

(T-term in quenched lattice QCD with Wilson quark action on a 12 3 x 20 
lattice at (5 = 5.7. The latter is evaluated with the aid of the variant wall 
source method, which was previously applied successfully for extraction of 
tt-tt scattering lengths and rj' meson mass. We found that the disconnected 
contribution is about twice larger than the connected one. The value for 
the full ir-N a term a = 40 — 60 MeV is consistent with the experimental 
estimates. The nucleon matrix element of the strange quark density ss is 
fairly large in our result. 
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Study of the tt-N a-term defined as the nucleon matrix element of the scalar 
density, 

a = mu + md < N \uu + dd\N >, (1) 

has a long history and still shows some controversy. Current algebra and PCAC 
relates the a term to the quantity S = f%D + (2m 2 ) where D + (2m 2 ) is the cross- 
ing even tt-N scattering amplitude at the Cheng-Dashen point t = 2/i 2 , and a 
dispersion analysis of tt-N scattering leads to the value £ = 64 (8) MeV [l[] MeV. 
On the other hand, treating the flavor symmetry breaking to first order, one finds 
a = cr /(l - y) with ct « 25 MeV and y = 2 < N\ss\N > / < N\uu + dd\N >@. 
The two estimates cannot be reconciled unless the nucleon has an unexpectedly 
large strangeness content y ~ 0.6. A possible resolution is that the variation of the 
tt-N amplitude from the Cheng-Dashen point t = 2fi 2 to t = 0, where the a term 
(|I|) is defined, is substantial, resulting in a smaller value a ~ 45MeV||. Combined 
with the suggestion]^] that one-loop chiral perturbative corrections raises the value 
of o"o to o"o = 35MeV, this would imply a more reasonable value y ~ 0.2 of the 
nucleon strangeness content. 

In principle lattice QCD provides a more direct and more satisfactory route 
for resolving the issue through a numerical calculation of the a term matrix el- 
ement. A serious technical obstacle, however, has been that evaluation of the 
disconnected amplitude of quark loops and nucleon propagators projected onto 
the zero momentum state shown in Fig. 1(a) would require a prohibitively large 
number of quark matrix inversions if one were to carry it out with the conventional 
method of point source. In this article we show that the variant of the method 
of wall source without gauge fixing, which has been successfully applied for cal- 
culation of 7T-7T scattering lengths|| and rf meson mass[[|], also works well for the 
disconnected contribution to the a term. Employing the Wilson quark action, we 
find an encouraging result, a = 43(7)MeV, from an estimate of the matrix element 
and the up and down quark masses, or a = 58(ll)MeV from that of m^a '/m 2 , 
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albeit within quenched QCD at a fairly strong coupling of (3 — 5.7 with the lattice 
spacing a « 0.14fm. 

Let us note that previous estimates of the connected contribution in quenched 
QCD yielded values in the range 15 — 25 MeV@. The derivative m q dm^ / ' dm q in 
quenched QCD, equal to the a term via the Hellmann-Feynman theorem, also gave 
similar values 0]. These results were considered reasonable; sea quarks are absent 
in the quenched approximation, and hence one expects a connec t e d ~ o"o = 25MeV. 
This already suggested that the disconnected contribution is substantial. Indeed 
we find (J disconnected connected = 2.23(52). Previous attempts for estimating the 
disconnected contribution in full QCD, made either indirectly through a compar- 
ison of the slope m q dm^ / dm q and the connected contribution || or through the 
noisy estimator for quark loopflHj, indicated that the disconnected contribution 
is as large as the connected one. However, both results suffered from significant 
errors so that reaching a reliable conclusion was difficult. We should add that the 
large disconnected contribution poses a problem regarding the strangeness content, 
which we shall comment upon toward the end of this article. 

To extract the nucleon matrix element of the scalar density S = uu + dd we 
calculate the ratio of the three-point function of the nucleon and scalar density to 
the nucleon two-point function, each projected onto the zero momentum state 0, 

R (t) = — —. -, . — > const. + Z q < N\uu + dd\N > t, 2 

v ; < N(t)N(0) > 

with Zs the lattice renormalization factor for the scalar density. The t = time 
slice is excluded from the sum in the numerator since we fix the t — time 
slice of gauge configurations to the Coulomb gauge to enhance nucleon signals. 
This procedure affects only the constant term in @. The connected amplitude 



(Fig. 1(b)) can be calculated by the conventional source method [ 11 1. To handle 
the disconnected piece (Fig. 1(a)) we prepare a quark propagator evaluated with 
unit source at every space-time site ( except for the t — time slice ) without 
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gauge fixing @; 

G(n / ,0= E G(n',t>",t"). (3) 

(n",i'VO) 

The product of the nucleon propagator and Z)(n',t'^o) Tr[C(n', t')] equals the dis- 
connected amplitude up to gauge-variant non-local terms which cancel out in the 
average over gauge configurations. The superior feature of this method is that it 
requires only two quark matrix inversions for each gauge configurations. We note 
that with the exclusion of the t" = term in ([3]), the Fierz mixing of the quark 
propagator (|3|) for the disconnected loop and nucleon valence quark propagators 
is automatically avoided. 

We have applied the method for the Wilson quark action in quenched QCD 
at (3 = 5.7 on a 12 3 x 20 lattice. We analyzed 300 configurations for the hopping 
parameter K = 0.160, 0.164 and 400 configurations for 0.1665, generated with the 
single plaquette action separated by 1000 pseudo-heat bath sweeps. In order to 
avoid contaminations from the negative-parity partner of the nucleon propagating 
backward in time we employ the Dirichlet boundary condition in the temporal 
direction for quark propagators. We used the relativistic nucleon operator iV = 
(*gC _1 75g)g. Errors are estimated by the single elimination jackknife procedure. 

In Fig. 2 the ratio R(t) is plotted for the connected and disconnected amplitudes 
as a function of t for the case of K — 0.160 and 0.1665. For the connected 
amplitude we observe a clear linear increase with very small errors of less than 
1-2% up to t ~ 10 as has been known from previous workp], |9|]. Signals are also 
reasonable, albeit worse, for the disconnected amplitude, and are consistent with 
a linear behavior in t up to t w 10. The slope of the disconnected amplitude is 
even larger than that of the connected one, indicating a substantial contribution 
from the disconnected amplitude to the a term. To extract the scalar density 
matrix element we fit the data to the linear form (0) with the fitting range chosen 
to be 4 < t < 9. In Table 1 we tabulate the fitted values of the matrix elements 
corrected by the tadpole-improved renormalization factor in the MS scheme at 
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the scale \i = 1/a given byfl2|l 

- log aa - 0.0098^1 \(\- %—\ 



Z s = (^1 + a s (^- log a/i - 0.0098J J [1-—), (4) 

where we used ajjg{l/a) = 0.2207 for a s . 

We present the quark mass dependence of the matrix element < N\uu+dd\N > 
in Fig. 3 as functions of the bare quark mass m q = (1/K — l/K c )/2 using K c = 



0.1694[T3|]. The disconnected contribution (circles) is about twice as large as the 
connected contribution and increases almost linearly toward small quark masses, 
while the connected contribution remains constant. Extrapolating linearly in the 
quark mass, we find < N\uu + dd\N >= 5.8(1.4) for the disconnected contribution, 
2.62(6) for the connected piece, and 8.6(1.4) for the sum at m q = 0. The value for 
the connected contribution is reasonably consistent with the value of the derivative 
dm^jdmq = 2.97(11) obtained by fits of the nucleon mass in the range K = 
0.160 — 0.1665 where our calculations are made. 

An estimate of the a term in physical units requires the value of the bare 



mass 



13] rh — (m u + md)/2 for the up and down quark and the physical scale of 



lattice spacing. From fits of the meson spectrum obtained on our set of gauge 
configurations, we find m\ = 2.710(32)m, m p = 0.5294(86) + 1.593(61)m, from 
which we obtain m = 0.00342(12) and a' 1 = 1.45(2)GeV using m w = 140MeV 
and m p = 770MeV. These values yield a = 43(7)MeV, showing an encouraging 
agreement with the dispersion theoretic estimate of the a term. Alternatively one 
can form the dimensionless ratio m^a/m^ at each quark mass and extrapolate 
to the chiral limit. The ratio turned out to depend little on m g , and a linear 
extrapolation gives m^ajm^ = 2.80(53) or o = 58(ll)MeV with the physical 
nucleon and pions masses. The difference between the two values of a originates 
from the fact that m^jm p = 1.48(4) at f3 = 5.7, taken for the present simulation, 
is larger than the experimental value 1.22. 

Let us turn our attention to the nucleon matrix element of the strange quark 
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density < iV|s,s|iV >. This matrix element receives contribution only from the 
disconnected amplitude. For an extraction of its physical value one needs to ex- 
trapolate the valence quark mass to the chiral limit keeping the strange quark mass 
fixed. In Table 2 we summarize the results for the matrix element < iV|,ss|iV > 
corrected by the tadpole- improved Z factor relevant for such an analysis. Re- 
sults of a linear extrapolation to the chiral limit for the valence quark are also 
listed in Table 2. Fitting the extrapolated values for three values of K s we obtain 

< N\ss\N > Kval =K = 2.93(93) - l.l(7.9)m s . 

The strange quark mass may be estimated by generalizing the relation = 
2.710(32)m to m 2 K = 2.710(32) (m + m s )/2 and using the experimental ratio 
m K/ m p — 0.64, which yields m s = 0.0831(30) or K s = 0.1648. We then obtain 

< iV|s,s|iV >= 2.84(44). Two quantities of physical interest that can be esti- 
mated from this result are the strange quark contribution to the nucleon mass 
m s < N\ss\N > jm N = 0.302(48) and the K-N a term a KN = (m + m s ) 

< N\uu + dd + 2ss\N > /4 = 0.310(37) = 451(54)MeV. While these results appear 
quite reasonable, we need to note that the value of the matrix element < iV|ss|iV > 
itself is large. Combined with our result < N\uu + dd\N >= 8.6(1.4) for the up 
and down quark, we find y = 2 < iV|,ss|iV > / < N\uu + dd\N >= 0.66(15), while 
a phenomenological estimate gives y ~ 0.20. We have not identified an apparent 
origins of systematic errors that lead to the result, perhaps except for a possible 
scaling violation as our simulation is made at a rather strong coupling. This point 
should be examined through a repetition of the calculation with a smaller lattice 
spacing, which we leave for future investigations. 

To summarize we have shown that the method of wall source without gauge 
fixing 0, U helps to overcome the computational difficulty of yet another quantity, 
the disconnected contribution to the tt-N o term, for which we obtained o = 
40 — 60MeV with a 15% error at a lattice spacing a « 0.14fm. Our results also 
show a large strangeness content in the nucleon. 
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Figure Captions 

Fig. 1 (a) Disconnected and (b) connected contributions to the nucleon-scalar den- 
sity three-point function. 

Fig. 2 Ratio R(t) for the disconnected (circles) and connected (triangles) ampli- 
tudes at P — 5.7 on a 12 3 x 20 lattice. Solid lines are linear fits over 4 < t < 9. 
(a) K = 0.160 and (b) K = 0.1665. 

Fig. 3 Disconnected (circles) and connected (triangles) contribution to the matrix 
element < N\uu + dd\N > as a function of the quark mass m q = (1/K — 
l/K c )/2 in lattice units. Solid lines are linear fits for extrapolation to the 
chiral limit m q = 0. 
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Tables 



K 


0.160 


0.164 


0.1665 


K c = 0.1694 


m q 


0.173 


0.097 


0.051 




# conf. 


300 


300 


400 






0.85 


0.74 


0.60 




m N 


1.2957(79) 


1.080(10) 


0.926(13) 


0.783(16) 


< N\UU + dd\N > CO nn. 


2.323(15) 


2.413(30) 


2.693(81) 


2.615(61) 


< N\uu + dd\N > disc . 


3.56(76) 


4.58(92) 


5.14(1.10) 


5.8(1.4) 



Table 1: Disconnected and connected contributions to the a term matrix element 
< N\uu + dd\N > at (3 — 5.7 on a 12 3 x 20 lattice in quenched QCD. Results are 
corrected by the tadpole-improved Z factor. Values at K c are obtained by a linear 
fit in 1/K. 





< N\ss 


\N > 






K s = 0.160 


K s = 0.164 


K s = 0.1665 


0.160 


1.78(38) 


2.02(38) 


2.20(39) 


0.164 


2.04(47) 


2.29(46) 


2.49(45) 


0.1665 


2.50(55) 


2.81(55) 


2.57(55) 


K c = 0.1694 


2.67(69) 


2.99(68) 


2.77(68) 



Table 2: < iV|ss|iV > as a function of the hopping parameters of valence and 
strange quarks obtained at /3 — 5.7 on a 12 3 x 20 lattice in quenched QCD with 
300 gauge configurations.. Results are corrected by the tadpole-improved Z factor. 



11 




0.0 0.1 0.2 



m =(l/K-l/K)/2 



Fig.3 



This figure "figl-l.png" is available in "png" format from: 
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